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The Effect of Stochastic Noise on Quantum State Transfer 
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FIG. 1: Schematic arrangement of cavity and three level 
system. The small circles represent Faraday isolators so the 
forward and reverse propagating paths are distinguishable. 
Adapted from Cirac et al. p|]. 



T.M. Stacef] and C.H.W. Barnes 
Cavendish Laboratory, University of Cambridge, Madingley Road, Cambridge CB3 OHE, United Kingdom 

(Dated: February 1, 2008) 

We consider the effect of classical stochastic noise on control laser pulses used in a scheme for 
transferring quantum information between atoms, or quantum dots, in separate optical cavities via 
an optical connection between cavities. We develop a master equation for the dynamics of the 
system subject to stochastic errors in the laser pulses, and use this to evaluate the sensitivity of the 
transfer process to stochastic pulse shape errors for a number of different pulse shapes. We show 
that under certain conditions, the sensitivity of the transfer to the noise depends on the pulse shape, 
and develop a method for determining a pulse shape that is minimally sensitive to specific errors. 
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I. INTRODUCTION 

The ability to transfer a quantum state from one sys- 
tem to another will be a useful tool in the future of quan- 
tum information, both for quantum communication and 
for building large quantum computers from smaller com- 
ponents. Various schemes have been proposed for trans- 
ferring quantum states [Jl[ |[ || g, 0] , some of which 
are based on teleportation and entanglement purification. 
Others are more direct, having a quantum channel (typ- 
ically an optical fiber) connecting the two systems, and 
classical control pulses applied to each system ensure that 
emission from the first system and absorption by the sec- 
ond are perfect. 

Over long distances, it appears that protocols based 
on teleportation and purification are necessarily more ef- 
ficient than direct schemes ||, but over short distances 
(e.g. for interconnects between quantum components), 
the more direct method, first proposed by Cirac et al. 

may be practical, and is probably easier to imple- 
ment since it does not require quantum measurement. 
The scheme by Cirac et al. Q implements quantum state 
transfer from one atom located inside an optical cavity 
to another atom in a separate but connected cavity, me- 
diated by Raman transitions via an atomic state (outside 
the computational basis) as well as cavity and external 
photons, and is controlled by laser pulses incident on each 
atom. The transfer is therefore subject to amplitude and 
phase errors in the control pulses. 

The idea j^] is also applicable to single-electron spin 
based systems in semiconductors, with additional uses 
such as spin measurement [|J and as a controllable 
source of decoherence for initial state preparation. The 
intermediate state required for spin-flip Raman transi- 
tions between single electron spin states is a state known 
as a trion || |j| . This state is a bound state of two elec- 
trons and a single hole, and couples to the single electron 
spin states via x- and y-polarised light. Physically, this 
corresponds to the optical excitation of an exciton which 
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may form a bound state with the electron, very much 
akin to an H _ complex. Trion states have been observed 
experimentally in GaAs jlO| O, and are relatively long 
lived, and so the scheme has the potential to be realised 
in semiconductor systems. 



In this paper we present an analysis of the effect of 
two possible types of stochastic error in control pulses, 
and consider a method for reducing the sensitivity of the 
pulse to such errors. Our results show that it is possible 
to design pulse shapes that reduce the sensitivity of the 
transfer process to certain sources of errors in the pulse. 
We give a short account of the applicability of this scheme 
to a semiconductor device. 



In section II we briefly review the proposal for quantum 
state transfer given by and give a new analytic form 
for a laser control pulse s hap e that implements quantum 
state transfer. In section [IlJ we formulate a master equa- 
tion for the transfer subject to stochastic errors in the 
control pulse, and use the master equation to compute 
the sensitivity of the transfer process to such errors. In 
section |[y] we present a procedure to optimise the pulse 
shape to reduce the sensitivity of the transfer to control 
pulse errors. We then finish with some conclusions. 
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II. QUANTUM STATE TRANSFER 
A. Review of proposal 

In Cirac et al.'s proposal for quantum state transfer, 
0, two 3-level systems (atoms, in the original proposal) 
are located inside separate optical cavities, as illustrated 
in Fig. [I], with the left atom starting in an arbitrary su- 
perposition of its lower levels, \ip) — c g \g) + c e |e), and 
the right atom in the state \g). The objective of quan- 
tum state transfer is then to transfer the state from 
atom 1 to atom 2 by applying suitable laser pulses to 
each atom. This may be expressed as 

|^)ib)2®|0)i|0)a-^|fl)i|^) 2 ®|0)i|0)a, (1) 

where the subscripts indicate the atom/cavity system. 
The term |0)j is cavity j's photon state in the number 
basis, and is the vacuum state at the start and end of the 
process. 

In the case that the detuning between the laser fre- 
quency and the atomic transition from \g) to |e), lul—luq, 
is larger than all other frequencies in the problem, one 
can write a Hamiltonian for each atom-cavity system, 
labeled j = 1,2, from which the upper level \r) is adia- 
batically eliminated, 



Hj(t)/h = 



n 3 {tf 



loo) is the effective coupling strength between state \g) 
and |e), and is proportional to Ej(t), the amplitude of 
the laser electric field, which we will use to control the 
transfer. We have also made the following definitions: a,j 
is cavity j's photon creation operator, g is the cavity- 
atom coupling strength, fi er is the dipole matrix ele- 
ment between |e) and |r), e is the laser polarisation, and 
4>j(t) is the phase of the laser. In deriving this equation, 
there is an assumption that the laser frequency is detuned 
from the cavity mode , u> c , nearest the transition, by an 
amount u>i — u c — g 2 /{ujl — ^o), which is the ac-Stark 
shift of the state \g)j in the electric field of the laser. 



A(lu l - too) 



1 33^ 



igj(t)(< 



in 



Following the original proposal by Cirac et al. j|| , the 
transfer is made to be unidirectional and is constrained 
by the requirement that the photodetector, indicated in 
Fig. [l], never registers a photon. The entire system is then 
an open quantum system, and is described equivalently 
by a quantum master equation, or a quantum trajectory 
punctuated by quantum jumps p|. Since it will be- 
come important for the later development of this paper, 
we summarise now some important points from Cirac et 
al, who adopt the quantum trajectories approach to the 
problem. The non-Hermitian effective Hamiltonian that 
describes the evolution of the open system in between 
jumps is H c s = Hi+H2 — iK(a\ai + a2a2 + 2a2(ii), where 
Ma H ) K 1S ^ 6 cav ^y leakage rate and the jump operator is 
J c — ai+a2. In the language of quantum trajectories, the 



(2) 



Here, Q,- (t) = 2fi er ■ eEj(t)/H is the classical Rabi fre- 
quency of the atom-laser system and g(t) = g£l(t)/(u>L — 



system state vector evolves under the Schrodinger equa- 
tion with the effective Hamiltonian H e g. This Hamilto- 
nian does not increase the total excitation number of the 
system, so we write the state vector as 



J 



m)) = c g \gg)\00) + c e { ai (t)e^%g)\00) + a 2 (t)e^%e)\00) + \gg) (^j£(|01) + |10» + ^(|01) - |10» 

(3) 



r 



The initial conditions and final configuration (for a suc- 
cessful transfer) are 

oc\{— oo) = 1, 4>i(— oo) = and tt2(+oo) = 1, 02(+oo) = 0. 

(4) 

we refer to the condition at t = +oo as the zero-jump 
condition, which holds as long as as c\ip(t)) — for all t, 
and this requires that 

Ps{t) = 0. (5) 

If we choose the phase of the laser to vary with time 
according to <j>j{t) — Vlj(t) 2 / (4(u>l — u>q)) then all in- 
dependence factorises out of the Schrodinger, resulting 



in the linearly independent evolution equations 

ax(t) = gi (t)/3 a {t)/V2, (6a) 
d 2 (t) = -sa(*)/8«(t)/V2, (6b) 
j3.(t) = gi(t)ai(t)/V2 + g 2 (t)a 2 {t)/V2 + Kp a {t){6c) 

Finally, the normalisation condition for the state vector 
requires that 



\a,(t)\ 2 + \Mt)\ 2 + \Ps(t)\ 2 + \p a (t)\' 



1. 



(7) 



We could also write an evolution equation for j3 a , how- 
ever this is the same equation as obtained by differenti- 
ating Eq. (0) with respect to time, and then substituting 
Eqs. (pa|) and (6b), so is not linearly independent. 
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There are thus only five independent equations, 
Eqs. (§) through (@), but six unknowns, ai,a 2 , [3 S , (3 a , gi 
and g 2 , and so there is no unique solution that satis- 
fies the requirements of quantum state transfer. There 
is a notable symmetry in the evolution equations and 
the initial and final states which is that if the cavity la- 
bels j = 1,2 are exchanged and time is reversed, i.e. 
1 <-» 2 and t — » —t, then they are recovered, so long 
as g 2 (t) = gx(—t), which is referred to as the symmet- 
ric pulse condition. Under this condition, Cirac et al. || 
present a method for generating pulse shapes that sat- 
isfy the requirements of quantum state transfer. This 
procedure is as follows: 

1. Select gi{t) for t > 0. 

2. Solve Eq. ( |6a] ) and a reduced form of Eq. (|6c|): 

$a{t) = -9i{t)0a(t)/V2 - Kgi(t)ai (t), subject to 
the initial conditions ai(0) = k/ ^/2(gi(0) 2 + k 2 ) 
and/3 Q (0) = -^1-201(0). 

3. Define , 9l (i) =-(V2 K /3 a (|t|)+.gi(|t|)ai(|t|))/a 2 (|t|) 
for t < 0. 

Then gx(t) is defined at all times, and as long as 
<7i(+oo) > 0, this pulse shape will satisfy the require- 
ments of quantum state transfer. Cirac et al. || present 
a particular solution to the equations, which they show 
accomplishes the transfer successfully. 

We make the point that if gi(t) — + as t — ► +oo 
the method doesn't guarantee a successful transfer pulse 
shape, but it does allow us to evaluate the success of a 
transfer based only on knowledge of gi(t) for t > 0, sub- 
ject to the presumption that g\{t) : t < is determined 
from <7i(t) : t > 0, according to the above procedure. 



B. Alternate solution method 

Using Eqs. (|oa| ) and (^) to replace g\(t) and g 2 it) in 
Eq. (|6c|), along with the zero jump condition, Eq. (|^), 
as well as the normalisation condition Eq. (Q) results in 
the following reduced non- linear equation for a±(t) and 
a 2 (t), 



d_ 

dt 



ai (t) 2 - a 2 (t) s 



+ K(l-ai(t) 2 -a 2 (<) 2 ) = 0. 



(8) 



Thus, another approach to finding a suitable pulse shape 
is to invent another equation depending on ai(t) and 
a 2 (t) (i.e. /(ai(i), a 2 (t)) = 0) consistent with Eqs. (Q) 
and (0) and then solve this equation and Eq. d§|) to find 
ai(t) and a 2 (t), from which we can calculate the other 
functions. 

We can motivate an example of this method by mak- 
ing gi(t) — g 2 (t), so that pulses on each atom are the 
same. We then divide Eq. (§a|) by Eq. @, to find 
di(i) = — a 2 (t), and taking account of the initial con- 
dition, Eq. (^), we find that a\(t) + a 2 (t) = 1, which is 




FIG. 2: Plot of the quantities given by Eqs. (10a) to (10c 



of the form f(ai(t),a 2 (t)) = 0, and is consistent with 
Eqs. (§ and (@). Thus we can use this as a second equa- 
tion, along with Eq. ^ to solve for a\{t) and a 2 (t). By 
substituting for a 2 (t) in Eq. (||), we find that a±(t) sat- 
isfies 



di(t) + 2/s(ai(t) - ai(t) 2 ) = 0. 



(9) 



There are two solutions to this equation, a.\(t) — (1 ± 
e 2Kt ) _1 , however one is unphysical since it diverges at t = 
(violating the normalisation condition). The physical 
solution is ct\(t) = (1 — tanh(«;t))/2, and from this we 
can determine all the other quantities to be 

a 2 (t) = ai(-t) = (l + tanh(Ati))/2, (10a) 
/3 a (t) = -sech(/«t))/V2, (10b) 

gi{t) — g%{t) = k sech(nt). (10c) 



Notably, this solution has the property that gi(t) — > 
as t — * cxd, and it also satisfies the symmetric pulse 
condition, 92 (i) = i)j so is a member of the class 
of solutions generated by the method of Cirac et al. [0 . 
We now have two solutions out of an infinite range of 
possible ones (i.e. the solution presented in ||, and the 
one above), which we will compare in later sections. 



III. STOCHASTIC NOISE IN LASER PULSE 

Having presented an alternate pulse shape to accom- 
plish quantum state transfer, we now use the two solu- 
tions, that in Cirac et al. [0 and in this paper, as test 
cases for evaluating the effect of stochastic noise in the 
laser pulses on the fidelity of the transfer. We reformu- 
late the problem in terms of a quantum master equation, 
as opposed to that of quantum trajectories j|, ||, 0, since 
we will eventually take averages over an ensemble of noise 
realisations, which must be done with respect to the den- 
sity matrix rather than the state vector. 
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FIG. 3: Structure of relati onsh ips between density matrix 
elements arising from Eq. (11a), with respect to the basis 
B, which is indicated on the left. Elements within each 
sub-matrix are mutually coupled. Arrows indicate the di- 
rection of coupling between the sub-matrices, so for instance 
p5,5 depends on P3,s, but not vice- versa, i.e. derivatives are 
placed at the arrow head, for example the equation for ps,z is 
P5,5 = ff2(t) (p4,5 + P5a) ~ 2k (p 3 ,5 + p 5 , 3 + p Sl5 ). The arrows 
therefore indicate the flow of information, from the states with 
a single excitation in the left cavity to single excitation in the 
right cavity, but not in reverse, as determined by the one-way 
requirement used to formulate the equations. The bold ele- 
ments are those that may be non-zero at t = —00, and are 
given by \ip){ip\. 



A. Master Equation 

Since the process of information transfer is one-way, we 
adopt the cascaded quantum system approach, described 
in the text by Gardiner and Zollcr |l2| , to derive a master 
equation. The quantum master equation for the cascaded 
system here is given by 



P : 



KL{p}, 



---- \H x {t)+H 2 {t),p 

£{p} =2ciipa\ — a\a\p 
+ 2a 2 pa\ 
- 2{[a\,aip] + [pa\,a 2 ]} 



pa\a\ 



a 2 a\p - pa 2 a\ 



(11a) 
(lib) 



where p is the system density matrix. Subsequently, we 
leave out explicit time dependence notation from the 
Hamiltonians, H\ and H 2 , unless required for clarity. 
In lieu of writing out all 25 coupled, first-order differen- 
tial equations arising from the master equation, we give 
a schematic representation of the coupling between ele- 
ments of the density matrix in Fig. ^, where arrows indi- 
cate the direction of the coupling. From this schematic 
representation we can see how the equations ensure that 
information does flow unidirectionally from the left cav- 
ity to the right. The initial condition for this master 
equation is simply 

p(-co) = (|^) 1 | ff ) 2 |00))( 1 <^| 2 (.g|<00|). (12) 

For future reference, if a matrix representation of an op- 
erator is given, it will be expressed with respect to the 
basis B = {\g0g0},\e0g0),\glg0) 7 \g0e0),\g0gl)}, where 



the state vector notation represents the state of the left 
atom, left photon number, right atom and right photon 
number in order. The solutions we are interested in will 
be those for which quantum state transfer has been suc- 
cessfully implemented, so 



p(+oo) = (|s)iMa|00»(i<ff| a M<00|). 



(13) 



We will use the notation po(t) to indicate any solution to 
the noiseless master equation, Eq. ( |lla| ), which satisfies 
both Eqs. ( p"2| ) and ([l3|). For an arbitrary time t, po(t) 
will depend on the pulse shape, nevertheless all permis- 
sible solutions will agree at t = ±00. Of course, for those 
pulses that do satisfy the requirement of quantum state 
transfer, the master equation will give the same result as 
solving the quantum trajectory equations, Eq. 



B. Adding Stochastic Noise 

We add the effect of stochastic noise in the pulse shapes 
to the dynamics of the problem, and so we follow a similar 
development to that presented in recent work |e| [lJJ. 
We assume that for each realisation of a transfer attempt, 
the actual classical laser pulse, (jj(t), applied to atom j 
is perturbed from the desired shape according to gj (t) = 
9j(t) + £j(*)A<?j(t), wnere is a stochastic function 

of time, and Agj(t) is some function which will depend 
on the origin of the stochastic error (i.e. noise), and 
is considered to be small. For instance, if the noise is 
associated with fluctuations in the amplitude of the pulse, 
then Agj(t) = gj(t), so that §j(t) = (1 + €j(t))gj{t), and 
if it is associated with errors in timing between the pulses 
applied to the two cavities then Agj(t) — gj(t), both of 
which we consider later. This modifies the atom-cavity 
Hamiltonian from that with an ideal pulse, which to first 
order in £j(t) is 



H j (t)->H j (t)+£ j (t)h j (Ag j (t)). 



(14) 



The operator hj is Hermitian and depends on Agj{t). 

With the modified Hamiltonian given in Eq. (|l4|), the 
master equation reads 



Hi +Zi(t)Ki{t) + H 2 + &(t)h 2 (t),p 



nL{p}, 
(15) 



which we formally integrate on both sides to find p(t) 
and substitute back into the right-hand side (RHS) of 
the equation, to generate the second order term in the 
Dyson series. We then take averages over an ensemble of 
noise realisations, and assuming the one-time correlation 
terms have zero mean, i.e. (£j(t)) = 0, we arrive at the 
result 
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(p(t)) = - i/n [Hi + H 2 , (p(t))} + K L{(p(t))} 







hi, 


f dt'{ 




1 — OO 






hi, 


f dt'{ 




J — OO 











) 


+ 


hi, 


f dt'{ 

I — OG 










) 


+ 


hi, 


f dt'{ 

1 — OO 



(16) 



We write the averaged solution to the noisy master equa- 
tion, Eq. @, as a correction to po(t), the density matrix 
in the absence of noise, 



(p(t)} = po(t) + eiApi(i) + e 2 Ap 2 (t), 



(17) 



where ej is a measure of the (small) variance of the noise 
on laser pulse j, and will appear as a coefficient in the 
two-time correlation function (£'))• ^ n the case 

that the stochastic terms on each pulse are uncorrelated 
with one another, (t;i(t)£i(t')) = 0, which we will here- 
after assume , we can substitute Eq. (|l7|) into Eq. jl^ ) 
and use Eq. (11a) to derive two independent master equa- 
tions for the first order (in ej) correction terms Api(t) 
and Ap 2 (t), 



A Pj 



i/h Hi + H 2 , Apj + KL{Apj} 



+ (i/h) 2 



dt' 



(18) 



The quantity tjApj (t) is therefore interpreted as a correc- 
tion to the noiseless density matrix due to stochastic er- 
rors in pulse j. We note that the homogeneous part of the 
above equation is the same as the original noiseless mas- 
ter equation, Eq. ( 11a ), but it now has an inhomogeneous 
driving term that depends on the solution to the noiseless 
equation, po(t). We are interested in errors in the den- 
sity matrix associated with stochastic errors in the pulse 
shape, so we assume that the initial state of the density 
matrix can be prepared exac tly, i.e Apj(— oo) = 0. The 



solution to equation Eq. (11a) as t — ► oo is determined by 
Eq. (|l3|), but the solution to Eq. ( |l8| ) is not constrained 
and may depend on both gj(t) and Agj(t). 

In numerical investigations of the error, we assume that 
the initial state of the left atom is \tp) = |e), so that the 
only non-zero matrix element is (e0g0\po(— oo)|e0g0) = 1. 
This assumption is not restrictive, since as may be seen 
from the form of the coupling of the evolution equations 
in Fig. ||, the only initial condition that influences the 
matrix elements Pi,j(t),i,j > 2 is that on the element 
P2.i(—co) — (e0gQ\po(—cc)\e0g0) = 1, and by the lin- 
earity of the equations we obtain the solution (up to a 
scaling factor) for any non-zero initial condition. 

At the end of a noisless transfer, the only non-zero el- 
ement will be ((70e0|po(+oo)|<?0e0) = 1, but at that time 
the correction to the density matrix may be non-zero, and 



the corresponding element rjj(t) — (g0e0\Apj(t)\g0e0) 
(which depends on Agj(t)) we call the noise sensitivity. 
This is motivated by Eq. (n7|) , since errors in pulse j will 
perturb the final result away from unity 



(g0e0|p(+oo)|.g0e0) = 1 + e^+oo). 



(19) 



We can conclude from this that rjj(+oc) < 0, since €j > 
and tr(p) = 1. The relation above makes evident the fact 
that the stochastic (classical) noise is a source of deco- 
herence for the transfer, since the probability of finding 
the system in the state (gOeOl^+oo^gOeO) is less than 
unity. 



C. White Amplitude Error 

We first apply the above result to the case where the 
stochastic errors arrive from unwanted fluctuations in the 
amplitude of the pulse, as mentioned earlier. In princi- 
ple this could arise from intrinsic laser amplitude noise, 
although this is typically negligibly small. More likely 
sources of error are in the details of the pulse shaping 
and might arise from phase and amplitude errors of the 
modes selected to generate the pulse. In any case, we do 
not consider these issues here, but instead assume that 
the error has the property that Agj (t) — gj (t) , so that 



(20) 



We will also assume that the stochastic terms are 8- 
correlated, 



(^m j (t'))=e j S(t-t'). 



(21) 



With these assumptions the master equation for the cor- 
rections to the density matrix reduce to 



A Pj 



Hi + H 2 ,A Pj 



H{A Pj } 



1 

27? 



h , 



h 3 ,po 
(22)" 



where the factor of 1 /2 in front of the double commutator 
arises from integrating only half of the Dirac-<5 function, 
and Eq. (|2|) is of Lindblad form. 

We have solved the equations numerically (using 
the numerical differential equation solver provided with 
Mathematica 4.0) for two different pulse shapes, one 
given in this paper, Eqs. (fLCM, and the other taken from 



FIG. 4: The evolution of the noise sensitivities, rjj, for am- 
plitude errors. The upper curves (solid') correspond to the 
pulse shapes given in this paper, in Eqs. (lid), whilst the lower 
curves (dotted) are for the pulse shapes presented elsewhere 
M. The darker traces correspond to the solution for r]i(t) 
lighter traces are for 772 (t) . 



the proposal for quantum state transfer by Cirac et al. 
|^| . The results for the noise sensitivities are plotted in 
Fig. [|. For each pulse shape the noise sensitivities 771 (t) 
and 772(0 appear to converge to the same value as t — > 00, 
however they do differ by a small but definite amount be- 
tween pulses shapes, being —1 (to the numerical precision 
of the equation solver) and approximately —1.17. 



D. White Timing Error 

Using the same methods, we can consider the effect of 
timing errors between the pulses applied to the cavities. 
In this case it makes sense to define only one noise sensi- 
tivity, 77, since we may assume that the pulse on the first 
atom defines the origin of time, and thus it will be errors 
in the timing of the second pulse relative to the first that 
will contribute to state transfer errors. 

For timing errors, the noisy pulse shape is given by 



FIG. 5: The evolution of the noise sensitivities, r\, for timing 
errors. The lower curve corresponds to the pulse shapes given 
in this paper, in Eqs. ([jl]), whilst the upper curve is for the 
pulse shapes presented elsewhere pi. 



IV. OPTIMISING AGAINST NOISE 

From the two examples presented above, we can see 
that different pulse shapes are more or less sensitive to 
errors in their generation, depending on the nature of the 
error. We will now examine how this variation may be 
used to optimise the pulse shape against sources of error. 

The utility of this is in recognising that stochastic noise 
in control signals contributes to decoherence, which has 
been noted previously [Q. Whilst in that work, the 
decoherence attributable to stochastic noise was found 
to be much less than that due to other effects (environ- 
ment), the calculation was based on a nominally constant 
control signal with stochastic noise superposed. The sit- 
uation we consider here, though in a different system, is 
slightly more general in that the Hamiltonian is explic- 
itly time dependent. Thus the errors in the control pulse 
will not be simply the intrinsic (e.g. thermal) noise as- 
sociated with the power source, but also dynamic errors 
associate with signal generation, e.g. due to digitisation 
of a nominally continuous pulse, or errors in timing. 



&(*) = S2(t + £(*))> 

* <»(*)+ £(*)3a(*)> (23) 

where the second line follows by expanding to first order 
in the quantity It may be seen from this form that 
^9z(t) — S^COj an( i f° r the pusposes of illustration, we 
will again take the errors to be ^-correlated, = 
eS(t-t'). 

The result of the numerical solutions for the corrections 
to the density matrix are shown in Fig. ||. In this case 
the pulse shape presented in this paper is very slightly 
more susceptible to timing errors than the pulse shape 
given by Cirac et al. [pj. 



A. Optimisation Strategy 

In principle, the problem of optimising the pulse shape 
against noise may be understood in terms of constrained 
optimisation. We wish to choose a pulse shape that min- 
imises the noise sensitivity, 77^ (+00), subject to the con- 
straint that it satisfies the condition of pulse transfer, 
i.e. 02 (+00) = 1. Conceptually, this would provide 
us with another equation to supplement Eqs. (||) which 
then makes the number of equations and unknowns equal, 
thereby specifying the pulse shape uniquely. In practise, 
we could formulate this problem in terms of Lagrange 
multipliers ]lq] , and derive the following constrained op- 
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timisation equations 
T>r)j(+oo) 



Va 2 {+oo) 
a 2 (+oo) = 1, 



0. 



(24) 
(25) 



where the operator T>/T>gi(r) represent the functional 
derivative with respect to the function gi(r), and A is 
an unknown Lagrange multiplier. We have introduced 
two more equations, and one more unknown, to bring 
the total number of each to seven. Since we have chosen 
.92(t) = gi(—r) we do not consider functional derivatives 
with respect to gi (t). Although we don't prove it for- 
mally, we find that rji (+00) = 772 (+°°) hi all numerical 
calulations we have performed, so we only consider vari- 
ational derivatives with respect to 771 (+00). 

In fact, computing either of the above functional 
derivatives requires the solution of a set of differen- 
tial equa tions that are of essentially the same form as 
Eq. (11a), since they are formulated by taking variational 
derivatives of the master equation with respect to g\ (r) , 
which must then be integrated. Since we do not have 
an analytic solution to the master equations for an arbi- 
trary pulse shape, it is difficult to see how practically to 
make use of these equations, except for very small scale 
discretisations. 

Another approach is to use numerical substitution. 
Here, we discretise the pulse shape, defining gj = gi(tj) 



for several discrete times 



tn 



! "ft 



(setting g n 



0), then interpolating between points. We allow all but g 
to be variational parameters, and go is chosen depending 
on gi, . . . ,g n -i to satisfy the constraint, «2(+oo) = 1, 
which is the numerical equivalent of a substitution for 
go in terms of all other variational parameters. We then 
have n—1 variational parameters, g\, ■ ■ ■ , g n —ii available 
for unconstrained optimisation of 771 (+00), for which we 
may use an unconstrained optimisation scheme. In the 
details of this method, the discretisation points are cho- 
sen to be in the domain t > 0, and go is determined from 
gi, ■ ■ ■ ,g n -i by use of the method outlined by Cirac et 
al. H (and at the end of section II A in this paper) which 
determines the efficacy of a given pulse shape in terms of 
its form for t > 0. 



B. Numerical Results 

A number of numerical optimisations are performed 
according to the numerical substitution scheme outlined 
above. In one set of optimisations, three points are al- 
lowed vary, and one of them, at t — 0, is used to sat- 
isfy the condition for pulse transfer. The locations of 
the other two points are equally spaced between t — 
and a fourth location, t — T, where the function is zero, 
i.e. the two intermediate points are located at t = T/3 
and t — 2T/3. These intermediate points are varied to 
optimise the pulse shape. The pulse shapes that result 
from the three-point optimisation process are presented 



as the lighter traces in Fig. |3(a) for five different values 
of kT = 2,4,6,8 and 10. 

In the other set of optimisations, six points are allowed 
to vary, again with the one at t = being used to satisfy 
the condition for successful pulse transfer. The remaining 
five points are also equally spaced between t — and 
t = T, for the same values of T as above, and the optimal 
pulse shapes for the six-point optimisations are shown as 



the darke r trac es in Fig. 6(a) 



Figure 6(b) demonstrates that the various pulses 
shown in Fig. 6(a) do each implement quantum state 
transfer successfully, since 1 0^2 (*) | 2 — * 1 as t — > 00 for 
all pulses. The figure also shows that as the width of 
the pulse becomes broader, proportional to T, the trans- 
fer process takes longer, so the slope of |a2(i)| 2 becomes 
shallower. 



Figure 3(c) show s the noise sensitivity, 771 (t) , for each 
Naturally, the value of 771 (+00) 



of the pulses in Fig. S(a 



gets closer to zero (i.e. better) as the number of optimi- 
sation points is increased, and it also gets better as the 
width of the pulse increases. It may be compared with 
Fig. |^ which shows the noise sensitivity for the two ana- 
lytic solutions discussed earlier and the optimisation pro- 
cess has clearly produced pulses that are significantly less 
sensitive to noise in the laser pulses than either of the an- 
alytic solutions discussed earlier in this paper. The best 
pulse shape computed in the optimisations performed 
here is that with T = 10, and has the limit 771 — ► —0.58 
as t — > 00. It appears that the optimisation procedure 
selects solutions with lower amplitudes, and this occurs 
because the form of noise we have assumed in this set 
of examples is proportional to the amplitude of the con- 
trol pulse, thus lowering the amplitude tends to lower the 
effect of the noise. 



Figure 6(d) shows that the photon wave-packet trans- 
mitted between the cavities, which is proportional to 
/3 a (t), becomes broader. It may be seen from Eqs. (0) 
and (U) that \ f3 a (t)\ 2 dt = 1, and the numerical solu- 
tions sati sfy t his to around one part in 10 5 . 

Figure |6(e)| plots the limiting value of 771 versus T, the 
time after which g\ is zero, for the three-point (light) and 
six-point (dark) optimisations. The solid curves are fitted 
using a function of the form a+ -^-pf, where a, b and c are 
fitting parameters. In both cases, the surprisingly close 
fit has a = —0.5 to better than one part in 10~ 3 , suggest- 
ing that the ultimate limit to the optimisation process is 
an infinitely wide pulse with 771(4-00) = —0.5. This indi- 
cates that the effect of the noise cannot be removed by 
allowing the pulse duration to become infinitely long. 



V. DISCUSSION 

The results presented in Fig. |^ above show that there is 
a clear possibility of selecting pulse shapes that are more 
or less sensitive to fluctuations therein. Figures I and I 
show that pulses which perform better with respect to 
one source of error may perform worse with respect to 
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2 4 /# S 10 12 



FIG. 6: (a) Optimised pulse shapes using three equally spaced points (light), and six equally spaced points (dark) for each 
of five different pulse widths, (b) Plot of a?2(i) 2 = (eO<70|po(i)|eO<70) for each of the pulses, (c) Plot of r/i(t) showing noise 
sensitivity for the same pulses, (d) j3 a {t) 2 = 2{glg0\po(t)\glg0) for each pulse, (e) Showing 771(4-00) versus T (points), and a 
fit to the data (solid) using a hyperpbolic relation. In all figures, traces from steepest to shallowest correspond to narrowest to 
broadest pulses in (a), with light traces corresponding to three-point optimisations, and dark traces to six-point optimisations. 



another, so in a realistic implementation of the optimi- 
sations presented here, detailed knowledge of the sources 
of stochastic error in the pulses would be required. 

The reduction in transfer error from an optimised pulse 
is moderate in the case of white amplitude noise consid- 
ered here, reducing the error in the transfer by a factor of 
a little over two compared to the original pulse shape pre- 
sented by Cirac et al. H . Though this is not a very large 
improvement, regardless of the shape chosen, pulse shap- 
ing will be needed to implement quantum state transfer, 
and so it is prudent to select the pulse with lowest noise 
sensitivity. 

The order of magnitude of transfer errors produced 
by the pulse are set by the magnitude of ej, which is 



essentially the relative variance of the pulse fluctuations, 
integrated over time. As mentioned earlier, for white 
amplitude noise, one candidate that could produce this 
is intrinsic amplitude fluctuations in the laser. For a 
cw laser, this can almost immediately be ruled out as a 
legitimate error source, since the relative RMS amplitude 
lasers in e.g. a diode laser are of order 10 -6 p6[ , so the 
variance is ej ~ 1CP 12 , which is completely negligible. 

A more likely source of error may be due to either the 
pulse generator or the pulse shaper. Whilst the orig- 
inal proposal for quantum state transfer was based on 
atomic states, we have had in mind the application of 
this scheme to electron spin states confined to quantum 
dots, with the dots contained within a linear microcav- 
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ity, in a similar scheme to that illustrated by Benson et 
al. and suggested by Imamoglu et al. || ||. For a 
cavity of length 1 /im, the free spectral range is around 
10 14 Hz, and if we imagine one of the Bragg mirrors has 
transmissivity of 10~ 3 , then the cavity loss rate is around 
k ps 10 11 sec _1 . This is an optimistic but not unreason- 
able estimate for the transmissivity, e.g. Lidzey et al. |l§| ] 
report a microcavity Q-factor of 125 (based on a cavity 
linewidth of ~ 20 meV and a mode frequency of 2.88 eV) 
for a A/2 cavity, which corresponds to a mirror trans- 
missivity of T m ; r = 0.025 |ll|, and good quality cavity 
linewidths may be as low as 0.5 meV [^0|, corresponding 
to T mir = 5x 10~ 4 . 

Since 1/k sets the time scale of the control pulse (see 
e.g. Fig, p), the typical width of the pulse is ~ 10 ps. 
Weiner ]21| gives a relation for the shortest temporal fea- 
ture that one may generate as St = T/tj where T is the 
largest temporal window, and r\ is a measure of the poten- 
tial complexity of the pulse. Typical values are quoted 
as T = 26.4 ps, which is the same order of magnitude 



as 1/k, and rj — 264, so St 
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-13 



s. Thus, if we as- 



sert that fluctuations at time scales less than St are un- 
controlled, then an upper estimate on relative amplitude 
fluctuations might be Sg w g'(t)6t, which means that un- 
controlled relative amplitude fluctuations from the pulse 
shaper could be as large as \ jr\ ~ 0.01, which is certainly 
not negligible. 

We note in passing that the time scale of 10 ps is 
roughly coincident with the characteristic time required 
for an electron in a travelling quantum dot generated by 
a surface acoustic wave (SAW) (velocity w 3000 m s _1 ), 
to traverese a distance of 1 /im, which is the typical di- 
mension of a SAW wavelength p2| , as well as a typical 
minimum spot size for an optical CW laser at its focal 
point. Therefore, if a microcavity were constructed at 
the end of the SAW based quantum computer proposed 
in , then a CW laser focused at the center of the mi- 
crocavity with a cross-sectional intensity profile shaped 
such that an electron passing through the microcavity ex- 
periences a time varying electric field which implements 
quantum state transfer (e.g Eq. (10c)), then the quan- 
tum state transfer scheme discussed in this paper could 
be used for both an interconnect between several SAW 
quantum processors and for electron-spin measurement 
(detection of a photon out of the cavity corresponds to a 



definite spin state of the electron). 

This paper has been fundamentally concerned with 
time dependent Hamiltonians, which are, conceptually, 
generally applicable to schemes for implementation of 
quantum computation. We therefore suppose that in 
many quantum gates implemented with time dependent 
control pulses, there could be some advantage in shaping 
the control pulses to be minimally sensitive to errors, so 
the ideas discussed in this paper may well be generalis- 
able to other quantum information applications. 

VI. SUMMARY 

We have analyzed the effect of stochastic noise on a 
control pulse that implements quantum state transfer 
between two atoms or quantum dots. We have found 
an alternate analytic solution for laser pulses that imple- 
ment quantum state transfer, and using the two analytic 
solutions we have available, we found that although dif- 
ferent control pulses are capable of implementing quan- 
tum state transfer, they may have different sensitivities 
to noise, depending both on their shape and the source 
of the noise. We then provided a method for calculating 
pulse shapes with lower noise sensitivities thereby opti- 
mising the pulse shape. This method was demonstrated 
to produce more optimal control pulses for the case of 
noise that is proportional to the control pulse amplitude 
with a white spectrum. A number of successively bet- 
ter pulse shapes were presented, and extrapolating the 
improvement in the noise sensitivity to its limit (i.e. for 
infinitely wide pulses) strongly suggests that the influ- 
ence of this particular form of noise cannot be cancelled 
completely. 
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